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We develop a subtractive renormalization scheme to evaluate the P-wave NN scattering phase 
shifts using chiral effective theory potentials. This allows us to consider arbitrarily high cutoffs 
in the Lippmann-Schwinger equation (LSE). We employ NN potentials computed up to next-to- 
next-to-leading order (NNLO) in chiral effective theory, using both dimensional regularization and 
spectral-function regularization. Our results obtained from the subtracted P-wave LSE show that 
renormalization of the NNLO potential can be achieved by using the generalized NN scattering 
lengths as input — an alternative to fitting the constant that multiplies the P-wave contact interaction 
in the chiral effective theory NN force. However, in order to obtain a reasonable fit to the NN data 
at NNLO the generalized scattering lengths must be varied away from the values extracted from the 
Oh' so-called high-precision potentials. We investigate how the generalized scattering lengths extracted 

from NN data using various chiral potentials vary with the cutoff in the LSE. The cutoff-dependence 
of these observables, as well as of the phase shifts at Tiab ~ 100 MeV, suggests that for a chiral 
■ potential computed with dimensional regularization the highest LSE cutoff it is sensible to adopt 

is approximately 1 GeV. Using spectral-function regularization to compute the two-pion-exchange 
potentials postpones the onset of cutoff dependence in these quantities, but does not remove it. 

PACS numbers: 12.39.Fe, 25.30.Bf, 21.45.-v 
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', theory is based on the spontaneously broken chiral symmetry of QCD. In it the long-range part of the strong interaction 
between nucleons is given by the exchange of QCD's approximately massless Goldstone bosons — the pions. This 
_ ^ dynamics is described by an effective Lagrangian which inherits the SU{2)l x SU{2)fi symmetry of QCD — as well as 
I the pattern of its breaking — and so reproduces QCD in the energy region below the chiral-symmctry-brcaking scale 
' A^sB ~ 1 GeV. In order for the effective theory to accurately represent QCD the low-cncrgy constants that encode 
the effects of degrees of freedom with energies > A^sb must be determined from cither lattice QCD calculations or 
experimental data. 

Since this theory is a quantum field theory, it requires renormalization. But, as an effective theory, it contains 
5J] ■ non-renormalizable operators. As a result one needs to establish a power-counting scheme in order to renormalize 
\ the theory order by order. The usual chiral counting in the single-nucleon sector involves a perturbative expansion in 
powers of Q ( "chiral perturbation theory" ) where Q is cither the momentum of the particles involved or the pion mass: 
Q = {p, m7r}. Weinberg [H, showed that a similar x^T expansion could be made for the NN potential, although 
infra-red enhancements due to the presence of two heavy particles rendered such a perturbation theory questionable 
for the NN amplitude. Since then considerable efforts have been devoted to the construction of NN potentials that 
are based on a XPT expansion [1, i, i, i, , with calculations now complete up to N'^LO {=0{Q'^)). These potentials 
have both long-distance [r ~ 1/m^) and short-distance (r <^ l/rriTr) parts. A key feature of the resulting "chiral 
effective theory" is that only certain short-distance structures are permitted at a given order in the Q expansion, 
because a momentum expansion is being employed. 

With an NN potential available it would seem straightforward to iterate that potential in a Lippmann-Schwinger 
equation (LSE) and so obtain the NN scattering amplitude, and from there phase shifts [sl. [sl. IgI. YA IsI. [qI. IToI. [ill. . 
How to perform such a non-pcrturbative treatment consistently remains an open question. It requires a cutoff 
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I. INTRODUCTION 

In the present work we will focus on nuclcon-nuclcon (NN) scattering in chiral effective theory (xET). This effective 
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in the LSE, because the xPT potentials grow with the exchange momentum q. The theory cannot be regarded 
as "properly renormalized" unless it produces results that are cutoff independent — or contains only weak cutoff 
dependence that can be removed by including higher-order short-distance operators. A contrast to such use of 
potentials in the LSE is provided by attempts to establish a perturbative expansion for nuclear forces that incorporates 
chiral symmetry. But st raig htforward implementation of this perturbation theory results in poor convergence in the 
key ^Si-'^Di channel [l3l.ll4l. [T5l.[l6l| (c.f. the recent work of Ref. This failure of perturbation theory is associated 

with the appearance of a new scale in the NN problem — one associated with the tensor part of the NN force [13, [3 ■ 
Rcnormalization-group methods can then be used to organize and implement a mixed treatment [20l . [2ll . [2^ , where 
one-pion exchange is treated non-perturbatively, and all parts of the long-range force that are of higher chiral order 
are evaluated in distorted- wave Born approximation (see also Ref. H^). But, in any of these analyses where pion 
effects are treated non-perturbatively, the complexity of the potential means that all analyses of the residual cutoff 
dependence of results are necessarily numerical. 

The operators of lowest chiral order that appear in the NN potential are two S-wave contact interactions, together 
with the one-pion-exchange potential. These operators are all 0(Q°) and have been regarded as forming the "leading- 
order" potential. In Refs. H, it was shown the theory is properly renormalized in the ^So, ^Si— ^Di channels 
and Ref. [26i| found that the same was true in some of the higher partial waves ("'^Pi, ^^2, ^Fs, "'^G4, ^Pi, '^Fa) ^. In 
contrast, Ref. [1^ showed that xET is not properly renormalized at LO in spin- trip let channels where an attractive 
singular tensor force acts (e.g. '^Pq, '^D2, '^P2— '^F2) (see also Ref. [l^]). In Ref. [1^ it was suggested that a contact 
term needs to be added in these channels — even though this breaks the straightforward xPT counting for short- 
distance operators. This behavior in different NN partial waves can be understood as a result of the singular behavior 
of the xET NN potential, with the channels where the leading singularity of the potential as r ^ corresponds to 
an attractive force needing a counter term to stabilize the phase-shift predictions as a function of the cutoff [1^, [2^ . 
Partial waves where the leading singularity has a positive coefficient — and so the dominant effect at short distances 
is repulsion — have A- independent predictions at large A (29l . [30| . A renormalization-group analysis confirmed these 
findings [H]. Ref. [3l[ questioned the usefulness of conclusions based on a leading-order analysis, as well as the 
necessity of considering cutoffs > 1 GeV. 

In this paper we extend the analysis of the xET potential in the NN P-waves at large cutoffs to the NLO and NNLO 
potentials derived in Refs. [1, 0, @. The key issue we are concerned with is whether renormalization of the LSE can 
be successfully accomplished when these higher-order potentials arc employed. This problem has been addressed in 
Refs. j3(l, 32] using NLO, NNLO, and N'^LO potentials and a co-ordinate space analysis. However, Ref. [l^l focused 
on the ^So channel, and the analysis of Ref. [1^ only yields results for A — * cxo and cannot address the approach 
to that limit. We know from the detailed studies of Refs. [1, 0, 0| that there exists a range of cutoffs A = 600-800 
MeV where the renormalization of the integral equation can be successfully carried out. However, if this analysis 
cannot be extended to cutoffs that are large with respect to the chiral-symmetry-breaking scale we are forced to 
conclude that use of the xPT potential as a non-perturbative kernel in the Lippmann-Schwinger equation only makes 
sense at cutoffs A < Ajj-gB- For practical application of chiral potentials to problems in, e.g., nuclear structure or 
electromagnetic reactions, knowledge of the maximum sensible value of A in such a treatment is a crucial component 
of forming accurate error estimates at a given order. 

At next-to-leading (NLO) and ncxt-to-next-to-leading order (NNLO), the xPT potential acquires two-pion-cxchange 
(TPE) parts, as well as seven short-distance operators of chiral dimension two, four of which arc in one-to-one 
correspondence with contact interactions in the NN P waves Q . Here we will analyze the P-wave predictions due to 
these potentials using a method based on subtractive renormalization of the LSE [H, [H, Is^] . In this approach, we 
do not fit the constant that multiplies the short-distance operator in the NN potential directly to data. Instead, we 
perform manipulations on the LSE in order to obtain an integral equation where we can directly input a (generalized) 
scattering length an which is obtained from low-energy P-wave phase shifts. (An earlier, alternative subtractive- 
renormalization technique involves invoking the Born approximation, and consequently is unreliable for the highly- 
singular potentials being discussed here [36l.[37|.) Using this method we can predict P-wave phase shifts at cutoffs up 
to 10 GeV, and examine the cutoff dependence of these predictions. We then examine how far an must be varied in 
order to obtain a best fit to phase shifts below Tiab = 100 MeV for each cutoff. These numbers are compared with the 
ail's obtained from the low-energy phase-shift analysis, in order to see whether there is evidence for a critical value 
of the cutoff at which it is difficult to fit both the threshold data and the data in the range up to 100 MeV. We also 
perform the same calculation but replace the contribution (and later, the whole) to TPE by the spectral-function 
regularization (SFR) [H] results with a SFR cutoff A = 800 MeV therein. 



^ Here we do not discuss the issue of quark-mass dependence, although it should be noted that iteration of a potential of fixed xPT order 
also leads to inconsistencies as regards the counting of short-distance operators proportional to the quark-mass [11, [13. 
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The structure of our work is as follows. First, in Section |TT1 we introduce our subtractive renormalization scheme 
for P-waves. Then, in Section lllli we consider the results when o ne-p ion exchange alone is taken as the long-range 
potential in these partial waves, and reproduce the results of Ref. [20] with our subtractive renormalization method. 
In Section IIVI we list the TPE obtained via dimensional regularization (DR) and spectral-function regularization 
(SFR) and the relevant partial- wave decomposition formulas. In Section |Vl we use these potentials, together with the 
subtracted LSE, to obtain phase shifts for all four NN P waves over a range of cutoffs from 500 MeV to 10 GeV. We 
then discuss the implications of our results for xET in the NN system and summarize our findings in Section IVTl 



II. SUBTRACTIVE RENORMALIZATION IN P-WAVES 



In this section we will introduce an approach to the renormalization of the Lippmann-Schwinger equation that 
is based on subtractive renormalization |34 . [s^ . This has the advantage that only renormalized quantities appear 
in the final equation, with all reference to the coefhcient of the contact interaction, the NN low-energy constants 
(LECs), having disappeared. Such techniques are in general useful for potentials that are the sum of a long-range 
part and a contact interaction. The difference between Ref. [s^ and this work is that, in contrast to the constant 
contact interactions considered there, we now treat contact interactions that arc operative in P-wavcs, and thus are 
momentum dependent. 

The partial-wave LSE has the following form 

tMP,P;E)^vMP,P) + j:-Mj^ , (2.1) 

where Pq/M ^ E is the center-of-momentum (cm.) energy and A the cutoff parameter. The angular momenta of the 
incoming (outgoing) momenta are indicated by the indices I {I'). We consider potentials of the form: 

vl^,fip',p) = vllfip',p) + Cf,^p"p\ (2.2) 

where p(p') indicates the incoming (outgoing) momentum in the center-of-momentum (cm.) frame, and v^f' is the 
long-range potential that is operative in this channel. For the purposes of this work vl;^ consists of at least one-pion 
exchange, and, in most cases, includes two-pion exchange too. The second term is the simplest contact interaction 
that can appear after partial-wave decomposition in this channel. We explicitly express the dependence of this term 
on the quantum numbers S and J, indicating that for each spin and total-angular-momentum channel there is, in 
principle, a different constant. 

If operators up to 0{Q^) are retained in the xET potential then the four operators necessary to yield contact 
interactions in the NN P-waves (^Pi, "^Pq, '^Pi, '^P2— '^F2) arc present — and indeed are needed in order to rcnormalize 
the xPT loop diagrams in vflf- that diverge as q^. The piece of these interactions that survives after projection onto 
the I = I' = 1 partial wave is an operator proportional to p' • p, which produces the second term in Eq. (|2.2p in that 
case. 

The main idea of our subtraction method is to construct the fully off-shell partial-wave i-matrix from the knowledge 
of v^^ and the on-shell value of the t-matrix for zero energy. This procedure of subtracting LSEs at zero energy to 
cancel unknown constants was introduced in Ref. [s^ for s-waves and will be generalized in the following for partial 
waves with 1^0. 

First, a generalized scattering length for arbitrary angular momenta I and I' can be defined as 
Dividing the partial-wave LSE, Eq. (IIH), by p"- we obtain 



4J{p',p;E) ^ v§f{p',p) ^y2 M r'^ dp" p"^ vf,^„ip',p")tf,Jip",p;E) 
p'l'pl p'^'p^ ^-jf TT Jq Pq + ie — p"^ 



(2.4) 



Since vJ/P'(j)',p) ~ the ansatz of Eq. (|2.4p is general and can be applied to any partial wave. In the following 

we concentrate on P-waves (l = I' ~ 1). For ease of notation we drop for the derivation the indices referring to 
the angular momenta S and J, unless they are explicity needed. In the case of P-waves the full potential vfif {p',p) 
involves a contact interaction of the form C;?/ p'p. First we consider the half-shell and on-shell i-matriccs at zero 
energy: 
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Subtracting Eq. (j2.6p from Eq. (|2.5p and multiplying both sides by p', the unknown constant Ci'i cancels and we 
obtain 
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Here we have used that for P-waves il' = I = 1), hm^^o = i^- Note that for the coupled channels 

'^P2— '^F2, counter terms only exist for {l,V) = (1, 1), i.e., the '^P2— '^P2 part, if one considers operators up to 0{Q^). 
Thus, up to 0{Q^), the treatment for this channel is exactly the same as for the uncoupled channels, since there is no 
new unknown constant. All of the above limits are well defined and can be directly obtained from Eqs. (|4.ip - (|4.1ip 



in Section [TVl The only unknown left in Eq. (|2.7p is limj._fO 
techniques. 

With the half-shell t-matrix at zero energy, lim/j-^o 
obtain the fully off-shell t-matrix. 
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and one can solve the equation by standard 
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(2.9) 



Subtracting Eq. (|2.9|) from Eq. (|2.8p and multiplying both sides by p, the constant Cm again cancels and we arrive at 
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In the case I' = I we can then use the property of the t-matrix that tiii{k,p';0) = tiii{p' , k;0) to obtain 
linifc^oi *''' fc^'"^ ] and solve Eq. ([^4^ to find hldPjSlM, Qncc we obtain hldmM, or, equivalently, ti,i{p,p';0), 
we can proceed to calculate the on-shell t-matrix and the phase shifts using resolvent identities that connect the 
operator t{E) to the operator t{0). Those details are laid out in Refs. [ssl. [36l|. 



III. ONE-PION-EXCHANGE POTENTIAL IN P-WAVES 



In the chiral effective theory the operators of lowest chiral dimension which can contribute to the NN potential are 



the operators of the one-pion-exchange potential (OPE). The leading-order long-range potential is thus v 



LB, 
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which we consider in this section. It has the well known form 

where q = p — p'- When considering the OPE alone, we adopt the axial vector coupling constant gA = 1.29, 
and pion decay constant fj^ = 92.4 McV. (By choosing this value for gA in the OPE we are taking into account 
0{Q^) corrections to it that shift the ttNN coupling constant away from the Goldberger-Treiman result and give 
us the freedom to reproduce the empirical value.) In addition, the nucleon mass M ~ 938.926 MeV and pion mass 
TOtt = 138.03 MeV arc employed. In the S-wave channels, ^ Sp and -^Si— "^Di, the OPE can be properly renormalized by 
a constant contact interaction in each of the channels [I^, [2J, [2^ . These contact interactions represent the ignorance 
of contributions from high momenta, i.e. short-range physics. However, since the lowest order in which a contact 
interaction that is a polynomial of q can survive after partial wave decomposition is proportional to q'"'"' , there is no 
contact interaction in leading order in the P-waves unless one adds a higher-derivative operator to the leading-order 
potential, i.e. an operator that is 0{Q^) in the usual xPT counting in powers of q. 

Therefore, in the Weinberg prescription for computation of the xET potential, the OPE alone gives the LO phase 
shifts in P-waves. The vin in each partial wave can then be obtained using the projection formulae presented in the 
next section. The resulting phase shifts in the four NN P-waves are shown in the left panels of Figs.[l][2l Fig. [T] shows 
that the phase shift of ^Pi and ^Pi become independent of the choice of cutoff once A > 2000 MeV. However, for the 
attractive channels ^Pq and '^P2-^F2, the phase shifts in Fig. [2] are cutoff dependent. This reproduces the finding of 
Ref. [2^ , where the attractive nature of the tensor 1/r^ potential in these channels leads to strong cutoff dependence 
of the phase shifts. 

To stabilize these channels, a contact interaction must be included. Thus, we include contact interactions that 
survive in the P-waves after partial-wave decomposition. Once we do this, the subtractive renormalization described 
in the previous section can be employed and the LSE be solved. For our calculations we first employ the generalized 
I = 1 scattering lengths extracted in [4l| as the input value to our subtraction scheme. (See first two lines of Table [J) 
As shown in the right panels of Figs. [TJll these contact interactions absorb the cutoff dependence in the ^Pg and 
'^P2-^F2 channels. However, a strong cutoff dependence is now created in the ^Pi and '^Pi channels, and produces 
a resonance- like structure in the phase shifts. Note that since we adopt the generalized scattering lengths extracted 
in [U to perform the subtraction, the results are not necessarily the best fit to the Nijmegen phase shift analysis as 
we will discuss in more detail in Section fVl 

With the calculations shown in Figs. [T] and [2] we reproduce the finding of Ref. [2^: if cutoffs larger than w 600 MeV 
are to be considered then additional contact interactions must be included in the LO potential which act in the ^Pg 
and '^P2 channels. Adding those contact interactions in the other P-waves destabilizes the phase-shift prediction as a 
function of cutoff — as might have been expected given the fact that the OPE potential is either non-singular (^Pi) or 
singular and repulsive ('^Pi) as r ^ in these waves [29l. |30|. 



IV. TWO-PION-EXCHANGE POTENTIAL 



The two-pion-exchange potential (TPE), V2-k, appears first at ncxt-to- leading order (NLO) in the chiral effective 
theory. Here we will consider two parts of V2Tr'- the NLO (Q^) part and the next-to-next-to-leading order (NNLO) part, 
which is 0{Q^). In the NLO piece the divergent loop diagrams must be regulated in order to obtain the potential. 
This can be achieved by either dimensional regularization [1, 0, Q or spectral- function regularization [l^l- After 
regularization, counter terms up to 0{Q^) are added to the TPE to absorb the divergences. We first consider the 
TPE obtained via dimensional regularization (DR), and adopt the formulation of where the axial- vector coupling 
constant is gA = 1.26, and the pion decay constant /^r — 93 MeV. The ttN LECs that appear in the NNLO potential 
are given by ci = -0.81, C3 = -4.70, and C4 = 3.4 (ah in GeV"!) Q. We use the TPE as derived in Ref. [J], which 
includes the 1/M corrections. This is the correct form in the power counting M ~ 47r/7r A^^^-sb although a 
different counting is adopted in Refs. 0,113 ^• 



^ It could be argued that in order to achieve NNLO accuracy in this power counting we should also consider corrections to OPE. 

Ref. [39II points out that these corrections can be accounted for by adopting the potential of Ref. [J and modifying by an energy- 
dependent factor. Using the prescription of Ref. [39l | to include these corrections to vi^r in our calculation alters our phase shifts 
by only a small amount at the energies under consideration here. It has also been shown in Ref. [40l | that these l/M'^ pieces of wiir 
have only a small impact in the ^Sq channel. A precise description of phase-shift data is not our goal here, so we do not consider this 
correction further in what follows. 
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For the convenience of the reader the expKcit expressions for the TPE [J] are given below: 

V2tv = Vc + Ti • T2 Wc + \Vs +Ti ■T2 Ws] CTi ■ (T2 

+ [Vt + Ti ■T2 Wt] (Ti ■ q (T2 • q 

+ [Vso +Ti ■T2 Wso] i{cri + CT2) ■ (q X p) 

+ Wq + Ti • T2 Wq] (Ti ■ (q X p) (T2 • (q x p). 



(4.1) 



Here (Ti(ti) are the spin (isospin) matrices. The subscripts C, S, T, SO, and Q indicate the Central, Spin-Spin, 
Tensor, Spin-Orbit and Quadratic spin-orbit components of the potential. The parts that cannot be absorbed into 
contact interactions have the explicit form as listed in Ref. [1] . For the NLO V2Tr we have for the long-range (constant 
and terms excluded) part of the potential 
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and w = •\/4m^ -|- q"^ and q = |q|. Meanwhile, the 0{Q^) part of V2tt has contributions 
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^^o = %j^{m. + {Aml + q-)A{q)} 



(4.11) 



Note that although the NNLO potential behaves as |qp at large |q|. the counter terms that renormalize the TPE are 
only of 0(Q2). 
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Another way to regularize the infinite loop during the calculation of the TPE is to adopt a spectral-function 
regularization (SFR), which is designed to keep the contribution from the loop integral to energy-momenta that are 
in the region where xPT is manifestly applicable. The spectral-function regularized TPE has a form similar to that 
given above for the TPE obtained with DR, but the functions L{q) and A{q) in Eqs. (|4.2p - (|4.1ip are replaced by 

L^{q) and A^{q), which have the following form [38| : 



lHi) 

A^iq) 
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2q 



K^w^ + q^s^ + 2Aqws 



6 (A — 2m^)— arctan 
2q 
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(4.12) 
(4.13) 



where A denotes the SFR cutoff on the loop integrals that arise when evaluating the potential and s = y A^ — Am"^ . 

By defining the linear combinations Uk — Vk + (4/ — 3)Wk {K = C, S, T, SO, Q), with / being the total isospin, 
the partial- wave form of the potential can be written as Q: 



vi'i{p\p) = {l'sj\v{p',p)\lsj), 



(4.14) 



where l,l'{s) denote the angular-momentum (spin) quantum numbers. With this we list the different partial-wave 
potentials for the TPE d 
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V. RESULTS WITH DIFFERENT TWO-PION-EXCHANGE POTENTIALS 



In this Section we present the P-wave phase shifts obtained with the TPE of xET. They are obtained by solving the 
LSE via the subtractive renormalization introduced in Section |TT1 We organize our findings as follows: we will discuss 
three different versions of the TPE, each of which is supplemented by the OPE that was the focus of Section [1111 In 
Subsections A and B we consider the TPE obtained via dimensional regularization up to NLO and NNLO respectively. 
In Subsection C, we consider the NNLO TPE obtained by spectral-function regularization. For each of these three 
potentials we first consider the "bare" TPE, i.e. we solve the LSE with v — t^i^ +'y27r alone to see how strong the cutoff 
dependence of the phase shift is for this "un-renormalized" case. Then, we perform our subtractive renormalization, 
thereby incorporating contact interactions C;f/p'p'' into the potential. We first solve the problem by adopting the 
generalized scattering lengths af/ extracted in Ref. [ilj. Then we adjust af/ to find the best fit with respect to 
the Nijmegen phase-shift analysis [i^ . Throughout this work we only consider np phase shifts, since we are mainly 
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concerned with issues of renormalization of the NN potential obtained from xET. While the subtleties of isospin 
breaking can be computed systematically within this framework, they go beyond the scope of this work. 

The issue of whether a contact term is required for a singular potential to reach cutoff independence can be linked 
to the short distance behavior of the potential in the coordinate space (see, for example, Refs. [l^, [2^, H^). If the 
potential is singular and attractive at r — » 0, then the contact term is required, and if it is not then the phase shifts 
will have a stable A — > oo limit even in the absence of a contact term. In Subsection D, we summarize our findings 
for the different chiral potentials and link them to the potentials' coordinate-space behavior. 

A. Dimensionally regularized TPE at NLO 

According to the power counting in powers of Q, the TPE at NLO in P-waves should be associated with a contact 
interaction Cf/pp'. To obtain insight into its effect, we first adopt the "bare" TPE at NLO — i.e. we do not consider 
the contact interactions when solving the LS equation — and examine the A dependence of the resulting phase-shift 
predictions. For this case, the "bare" TPE contains the constant and terms excluded in Eqs. (|4.2p & (|4.3p . The 
importance of these terms diminishes as A is increased, with the effects becoming negligible once A > 3 GeV. 

The results are shown as a function of energy in Fig. [S] Plotting the phase shift in each channel as a function of the 
cutoff (see Fig.|31 black dotted line) makes it clear that the and '^Pi partial waves show strong cutoff dependence, 
whereas the ^Pi and "^Pq phase shifts converge for cutoffs A > 1200 McV. This feature is important because once a 
potential reaches cutoff independence in the LS equation adding further contact interactions may destroy this feature. 

Now, in order to explicitly see the effect of the counter terms Cf/pp', we adopt our subtractive method to solve 
the LSE. For this we need as input values the four generalized scattering lengths: limfc^o[^^^-^|p^] = ^fh- For these 
we first adopt numerical values that are between those of the sets given in the first two lines of Table lU That table 
lists the values obtained for the generalized P-wave scattering lengths from the Nijmll and Reid93 potentials [i^ . 
together with those for the Cd-Bonn [i^ and AVIS [i^ potentials. 

The results are shown in Fig. [5l and indicate that the phase shifts have less cutoff dependence than in the unrenor- 
malized case — especially in the Tiab < 50 MeV region. But it is not until A > 5000 MeV that most of the channels 
become cutoff independent, and the phase shift in the ^P2 partial wave does not even stabilize then. In particular, 
in the "'^Pi and '^Pq channels the contact interaction delays until A « 2000 MeV the convergence with respect to the 
cutoff that was obtained in its absence. Moreover, Fig. [5] makes clear that using the constants af/ extracted from 
the potentials of Ref. [ill together with the NLO xET potential does not yield a good description of the energy 
dependence of the Nijmegen phase shifts — even for laboratory energies Tiab < 100 MeV. 

In order to obtain a better description, we vary afi at each cutoff so as to match the Nijmegen phase shifts for 
Tiab < 100 MeV. This leads to the phase shifts shown in Fig. [51 In general, after performing our best fit we are able 
to match the Nijmegen P-wave phase shifts much better, especially for Tiab < 60 MeV. As we will show in the next 
subsection, the NNLO pieces of TPE are needed in order to improve the energy dependence of the predictions, and 
thereby allow a better fit at Tiab > 60 MeV. 

We now assess how our results change depending on whether aff is input from Table [T] or if it is adjusted to fit the 
Nijmegen phase shifts in the range Tiab = 0-100 MeV. If the results are significantly different then that means our 
effective theory has sizeable dependence on the choice of input — the renormalization-point dependence is large. In 
order to study this issue we denote the generalized scattering length which gives the best description of a phase shift 
as abest- We then plot Ubest against A in Fig. [7] and [51 For the case, the varation of abest with A is so large that 
it must be plotted separately. The huge deviation of abest{'^P2) from the Nijmll value of —0.28 fm'^ around A = 2500 
MeV is correlated with the lack of cutoff independence in that channel. The xET fails at this order when it is pushed 
to such cutoffs. 

But problems occur already at A w 1000 MeV. This issue of a difference between the abest necessary to fit phase 
shifts with the NLO TPE potential and the af/ obtained from "high-precision" potentials is particularly obvious 
in the ^Pi and ^Pq channels, where at A « 1400 MeV and A « 1000 MeV respectively, the phase-shift predictions 
diverge away from the Nijmegen phase shifts. At these particular cutoffs changes in the generalized scattering length 
cannot bring the xET phase shifts any closer to the "data" : the use of any finite aff as input results in almost the 
same phase shifts as a function of energy. (We have verified that the same results are obtained if a contact interaction 
Cfip'p is added to the potential and the coefficient is fitted to data.) 

This feature is correlated with bumps that are seen when the NLO results for the best-fit ^Pi and ''Po phase shifts 
at Tiab = 100 MeV are plotted versus cutoff (see black dotted line in the top two panels of Fig. [9]). Meanwhile, the 
lack of cutoff independence in the coupled ''P2-'^F2 channel is prominent in the fourth, fifth, and sixth panel of Fig. [S] 
The only wave here that has a stable phase shift at NLO after one subtraction with an adjusted to reproduce the 
low-energy "data" is the "^Pi. Over a cutoff range of A = 600-3100 MeV its 100 MeV phase shift varies by only a few 
degrees. 
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The above results suggest that the highest LSE cutoff one can adopt if one inserts the NLO TPE in the integral 
equation is 1000 MeV. As A goes beyond this value we encounter increasing difficulties in reproducing the Nijmegen 
P-wave phase shifts. 

B. The NNLO TPE computed with dimensional regularization 

We now adopt the "bare" DR TPE at NNLO and examine the resulting cutoff dependence of the P-wave phase 
shifts. The results are shown in Fig. 0] and Fig. [101 all of the P-waves exhibit strongly cutoff-dependent features 
when v^^ = vit^ + V2tt with V2-k computed up to 0{Q'^) using DR. Comparing this to the results obtained for the 
un-renormalized TPE at NLO in the previous subsection, we see that the NNLO part of TPE has caused the ^Pi and 
•^Po phase shifts to oscillate with respect to the cutoff too. 

Next, we add the contact interactions Cf/pp' to the potential. We again first adopt values of af/ based on the first 
two rows of Table m With this, the oscillation in the P-waves with respect to A is greatly reduced for A > 2000 MeV. 
However, as shown in Fig. 1111 most of the resulting P-wave phase shifts arc rather far away from the Nijmegen 
analysis. 

To remedy the situation, we again vary the four input values af/ to find the best fit for Tiab < 100 MeV and obtain 
the results shown in Fig. [T^l At the cutoffs presented in Fig. [T^l our best-fit results are generally comparable in 
quality to those presented in Refs. H, §• (Note that in 0], an exponential regulator is added to TPE, which further 
suppresses the TPE at momentum q w A.) Results for abest versus cutoff and the phase shifts at Tiab = 100 MeV for 
A = 600 - 3100 MeV are again plotted in Fig. [7] and Fig. respectively (red dashed line). For the ^Pi, ^Pq and ^Pi 
channels, there arc only minor variations of the phase shift with respect to A. Those variations correspond exactly to 
the variation of abest with A shown in Fig. [T] On the other hand, the ^P2 F2 phase shifts display oscillations with 
respect to A. At NNLO one subtraction is not enough to stabilize this channel. 

Figure [7] shows that changes of 20-40% from the values of af/ obtained from the 'high-precision' potentials are 
needed in order to achieve the best fit. In particular, abest shows oscillatory behavior for A = 800 — 3100 MeV in the 
•^Po, '^Pi and ^P2— ^F2 channels. The osciUations are « 10%, and appear to start or become larger once the cutoff 
increases above 1 GeV. Such changes in observable quantities with the cutoff are disturbing, but are not beyond the 
bounds of what one might expect due to N'^LO corrections. 

Now we examine the success of our fit of phase shifts as A changes. Fig. [T^ shows that there is no significant difference 
in the phase shifts for those chosen cutoffs within the region where the fit was performed, namely Tiab < 100 MeV. 
However, for Tiab > 100 MeV, the phase shifts computed with a cutoff A = 500 MeV differ from the ones with all 
other cutoffs for all four NN P-waves. This indicates that the cutoff of 500 MeV is too low: it cuts out too much of 
the TPE that enters the LSE. For A = 800 — 2000 MeV, the phase shifts show httle cutoff dependence apart from 
one specific case: the ^P2— ^F2 in the immediate vicinity of A = 1200 MeV. Here we are too close to the cutoff where 
these phase shifts diverge (see Fig. |9l). 

C. Spectral-function regularization TPE 

Both the phase shifts and abest should be (approximately) invariant over a range of cutoffs if we are to claim 
that the renormalization is done correctly. The oscillation of abest for A > A^ makes us re-think the role played by 
the NNL0((5^) part of the TPE in DR. This piece of the TPE contains terms proportional to jql"^, which are non- 
polynomial and have no corresponding counter term. Although the results in Fig. 1 121 show that the cutoff dependence 
of the phase shifts was absorbed by the 0{Q^) contact interaction, we still see that abest depends rather strongly on 
A in the "^Pq and "^Pi partial waves. Could it be possible that the reason for this is the existence of a piece of the 
potential that grows as |qp, and that docs not have a contact interaction to counter balance it? 

If we employ spectral-function regularization (SFR) to compute TPE and fix the cutoff on the spectral function, 
A, to a particular value, then the dominant power in the SFR TPE for |q| ^ A is |q|^. The long-range potential thus 
now has the same large-|q| as the contact interaction. Therefore, we will investigate how the phase shifts obtained 
when replacing the part of DR-TPE by SFR (with A = 800 MeV) diflfcr fr om those obtained in the previous 
Subsection. We focus particularly on whether a better convergence with respect to the LSE cutoff can be achieved. 
We call the resulting long-range potential in this case the "mixed" NNLO TPE. 

To see the effect of this "mixed" potential, in Fig. [13] we plot the P-wave phase shifts obtained by adjusting af/ 
to fit the phase shifts to the Nijmegen analysis for Tiab < 100 MeV. The most striking feature of the results is that a 
pattern similar to that observed in the DR TPE at NLO appears in the ^Pi and ^Pq channels, i.e., there is no way 
to describe the Nijmegen phase shifts for cutoffs near A = 1400 (1000) MeV for ^Pi(^Po). However, this is not very 
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surprising, since the leading singularity of these potentials is now given by the TPE in NLO, and thus the large-cutoff 
behavior is the same as for that potential. 

We do find that the oscillation in a^est with A that is present for the DR NNLO TPE does not appear when the 
"mixed" NNLO TPE is used: the variations of abest for all P-waves are under 5%. In addition, the values of abest 
are closer to the values obtained from the 'high-precision' potentials than are those foimd in the previous subsection. 
But the attempt to tame the part of the TPE at NNL0(Q3) by adopting the SFR NNLO pieces of the potential 
still fails for A > 1000 MeV, due to the issues in the ^Pi and '^Pq channels discussed in the previous paragraph. 

Finally, we adopt the full SFR TPE up to NNLO with the inner cutoff A = 800 MeV, and present the best-fit 
results in Fig. [9] (phase shifts at 100 MeV vs. A) as well as in Fig. [14] (phase shifts vs. energy at different cutoffs). 
The use of SFR removes the oscillation behavior in the phase shifts that occurs with the NNLO DR TPE. The phase 
shifts at Tiab = 100 MeV are now almost independent of cutoff from A = 600 to 2500 MeV. For A > 2500 MeV the 
"best-fitted" "^Pq phase shift at Tiab ~ 100 MeV becomes cutoff dependent. Thus, we conclude that for P-waves, once 
the SFR TPE is adopted up to NNLO, it is possible to perform the rcnormalization in the LS equation with a cutoff 
up to 2500 MeV. 

D. The co-ordinate-space connection 

The results of the previous three subsections are summarized in Table [ill In that table an entry "R" implies that 
a contact interaction is needed to make predictions in this channel independent of the cutoff. In contrast, an entry 
"U" implies that phase-shift predictions for that particular P-wave arc already A independent before any contact 
interaction is added. (See Fig. |4| for a graphical representation of this distinction.) If, in spite of this, a contact 
interaction is then added for A A^sb the phase-shifts either develop a resonant structure or become insensitive to 
the input value of af/ used to fit the coefficient of the contact interaction. The entries "U" correspond exactly to 
those cases where the leading r singularity of the potential has a positive coefficient, and so a repulsive singular 
potential dominates the physics at large cutoffs. Conversely, entries "R" occur in those channels where the leading 
r — > singularity of the co-ordinate-space potential corresponds to attraction [2^, [sO]- The ^P2 — ^F2 channel for DR 
NNLO is a special case, since the coupled-channels potential has two attractive singular eigenpotentials in the r ^ 
limit. Hence one subtraction is not sufficient to make phase shifts independent of A in this channel. The large- A 
behavior of our calculation can thus be completely understood from this co-ordinatc-space point of view. 

VI. SUMMARY AND CONCLUSIONS 

We have developed a subtractive rcnormalization scheme for xET NN potentials with contact interactions of the 
form Ciipp' and carried out calculations of the NN P-waves with one and two pion exchanges. We employed four 
different types of TPE potential: the dimensional-regularization (DR) TPE up to NLO and NNLO, the "mixed" 
potential (DR TPE up to NLO plus NNLO pieces based on spectral-function regularization), and spectral-function 
regularization (SFR) up to NNLO. Calculations using these xET NN potentials are done with and without the contact 
interaction. When the contact interaction is included our rcnormalization scheme only requires the knowledge of the 
general scattering lengths af,f , and thus allows us to go to an arbitrarily high cutoff in the LS equation. 

In the OPE case, we reproduced the results of Ref. [1^ with our subtractive rcnormalization scheme, i.e., for the 
attractive triplet channels (■^Pq and ^P2-'^F2) one needs to include a contact interaction that is nominally of higher 
chiral order so as to make the resulting phase shifts cutoff independent. Wc note that if such a counter term is added 
in all P-waves it destroys the cutoff independence in the ^Pi and ^Pi channels that is achieved with OPE alone. 

In case of the NLO (NNLO) TPE computed using DR we found that the generalized scattering lengths obtained 
from the Nijmll and Reid93 potentials in Ref. [4l[ produce phase shifts that are (approximately) cutoff independent 
after A > 5000(2000) MeV. The only exception is the '^P2-^F2 channel, where cutoff dependence persists to arbitrarily 
high cutoffs in both the NLO and NNLO cases. In the other channels, the final cutoff-independent phase shifts are 
rather far from those of the Nijmegcn analysis — especially at NLO. Wc therefore vary the generalized scattering to 
obtain the best fit for phase shift under 100 MeV. After doing that, we found that for the NLO TPE cutoffs above 
A = 1400(1000) MeV need to be excluded in the ^Pi(^Po) channel. At these cutoffs the NLO ^Pi and ^Pq phase 
shifts become insensitive to the input value of the constants afi- This occurs because the dominant piece of the 
co-ordinate-space potential as r — > in these channels is repulsive. Results at NLO also depend strongly on the 
rcnormalization point, since the energy dependence predicted there is rather unlike that of the data. 

At NNLO we are almost always able to fit all P-waves to the Nijmegen analysis — at least at the level of accuracy 
obtained in the NNLO fit of Ref. Q — after adjusting the values of the constants af/. In order to achieve such a 
fit we must take A > 600 MeV, as otherwise the LSE cutoff removes too much of the TPE physics that drives the 
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energy dependence of the phase shifts. But there is also an upper Hmit on the LSE cutoff at NNLO. There arc some 
cutoffs, e.g., A « 1250 and 2300 MeV, where the ■^P2 — '^F2 phase shift diverges. At this order one subtraction is not 
enough to render the ■^P2 — ■^F2 phase shifts stable with respect to cutoff. We would have to add another, e.g. tensor, 
contact term — a term that is nominally O(P^) — in order to produce cutoff-independent results. This is analogous to 
the situation in the ^Si - ^Di channel at NNLO and can be dealt with in the subtractivc approach provided an 
additional piece of experimental data is used as input [48l | . 

We then replaced the NNLO part of the DR TPE by SFR and again solved the LSE. Our results show that in this 
case the variation of atest with A is suppressed to less than 5%. A comparison between the results obtained with the 
NLO TPE and this "mixed" SFR TPE in the ^Pj-SPj channel implies that the mixed SFR NNLO TPE includes 
more physics than the TPE at NLO. Nevertheless, at larger cutoffs the repulsive NLO TPE potential in the "'^Pi and 
^Pq channels still leads to the same problem as observed in the DR NLO TPE, , and thus the allowed LSE cutoffs 
are stiU limited to A = 600-1000 MeV. 

Finally, we have the most success if we adopt the full SFR TPE at NNLO. In this case we obtain phase shifts which 
fit the Nijmegen "data" reasonablely well and are (almost) independent of renormalization point for LSE cutoffs as 
high as 2500 MeV. 

To summarize our study, let us define necessary conditions for renormalization of potentials in the LS equation. 
The resulting phase shifts should: 

(a) be (almost) independent of the cutoff A in the LSE, 

(b) fit the experimental analyses reasonably well, and be able to be improved order by order if the potential is obtained 
in a perturbative way, 

(c) be (almost) independent of the point of renormalization. 

Under these criteria, our best-fits with dimcnsionally-regularized chiral potentials arc limited to A < 1 GeV. Once 
cutoffs above this value are considered the Nijmegen phase shifts cannot be reproduced in certain channels at LO 
and NLO, while the instability of the '^P2 — '^F2 phase shift violates condition (c) at NNLO. Therefore we conclude 
that chiral potentials should not be inserted into the LSE if cutoffs A > A-^sb are employed. If spectral-function 
regularization is employed the criteria (a) and (c) are satisfied for A up to 2.5 GeV, but this is because another cutoff 
A, which itself is below A^^sb, has been used to ameliorate some of the features of the chiral potential that cause 
difficulties when it is inserted in the scattering equation. 

Note added: After this paper was submitted for publication we became aware of Rcf. [4^]. In that paper an 
analytically solvable model for NN scattering in S waves is presented. An EFT which has the same long-distance 
physics as the model is then constructed. Epelbaum and Gegelia show that when the EFT integral equation is solved 
with cutoffs larger than the EFT breakdown scale it leads to erroneous conclusions about the dependence of phase 
shifts on the parameters that determine the long-distance physics. Although the model of Ref. [40] does not contain 
the singular potentials discussed in this work, the practical implications for chiral EFTs of NN scattering arc in accord 
with the conclusions we have drawn here. 
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1.51 
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-2.41 


1.44 


-0.30 



TABLE I: The generalized P-wave scattering lengths gf/ as given in Ref. [4i[ for the Nijmll and Reid93 potentials, together 
with ones extracted from the Cd-Bonn and AV18 [4y| potentials, af-f is given in units fm^. 
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TABLE II; Need for renormalization in different NN P-waves when the xET potential is calculated to different orders. "U" 
implies that the potential does not need a subtraction to generate cutoff-independent phase shifts as A — > co while "R" means 
that a subtraction is required in that channel at that order. The * indicates that at NNLO the ''P2 — channel actually 
requires two subtraction to be made stable. 
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FIG. 1: (Color online) The NN phase shifts for the channels ^Pi and ^Pi as a function of the laboratory kinetic energy for 
different cutoffs A ranging from 0.8 to 5 GeV. In all cases the potential enters the LSE. The panels on the left 

show the results without subtraction, whereas the panels on the right show the results with one subtraction. The Nijmegen 
phase-shift analysis [ZJ is indicated by the open diamonds. 
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FIG. 2: (Color online) The NN phase phase shifts for the channels '^Po and '^P2— '^F2 as well as the mixing parameter £2 as a 
function of the laboratory kinetic energy for different cutoffs A ranging from 0.8 to 5 GeV. In all cases the potential 
enters the LSE. The panels on the left show the results without subtraction, whereas the panels on the right show the results 
with one subtraction. The Nijmegen phase-shift analysis [i?! ] is indicated by the open diamonds. 




FIG. 3: (Color online) The un-renormalized NN P-wave phase shifts as a function of the laboratory kinetic energy resulting 
from the dimensionally regularized TPE in NLO. The phase shifts are shown for cutoffs A ranging from 0.5 to 2 GeV. The 
Nijmegen phase-shift analysis [i^l is indicated by the open diamonds. 
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FIG. 4: (Color online) The un-renormalized NN P-wave phase shifts at Tiab = 10 MeV as a function of the cutoff for three 
different two-pion-exchange potentials: DR NLO, black dotted line; DR NNLO, red dashed line; SFR NNLO, solid green line. 
In each case v^^ = vin + vi-n ■ 




FIG. 5: (Color online) The NN P-wave phase shifts resulting from the use of OPE plus DR TPE at NLO with one subtraction 
as a function of the laboratory kinetic energy. Here the cutoff range shown is 0.5-10 GeV. The input value of aff is taken 
from Ref. [JJ] . The Nijmegen phase-shift analysis [i^l is indicated by the open diamonds. 
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FIG. 6: (Color online) The NN P-wave phase shifts as a function of the laboratory kinetic energy that result from choosing 
v'^^ = viTT + V2Tr, with the latter computed using DR TPE in NLO. Here the generalized scattering lengths afi are adjusted 
to give the best fit in the region Tiab < 100 MeV. The Nijmegen phase-shift analysis [13] is indicated by the open diamonds. 
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FIG. 7: (Color online) ajest versus A with the long-range potential chosen to be one-pion exchange plus DR NLO (black dotted 
line), DR NNLO (red dashed line) or SFR NNLO (green line). Here atest has been adjusted at each cutoff to give the best fit 
to the Nijmegen phase-shift analysis in the region Tiab < 100 MeV. The shaded region represents the range of the generalized 
scattering lengths extracted from the 'high-precision' potentials listed in Table I. 




FIG. 8: abeBt{'^P2) versus A for v^^ = vi^r +«2ir, with V2tt computed at NLO using DR. Here abest was adjusted at each cutofT 
to give the best fit to the Nijmegen '^P2 phase shift in the region Tiat < 100 MeV. 
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FIG. 9: (Color online) The renormalized NN P-wave phase shifts at Tiab ~ 100 MeV as a function of cutoff for three different 
V2n, with v^^ = vin + V2n. DR NLO (black dotted line), DR NNLO (red dashed line), SFR NNLO (green solid line). The 
generalized scattering lengths ctbest were adjusted at each cutoff to give the best fit in the region Tiab < 100 MeV (see Figs. [7] 
and HI). 




FIG. 10: (Color online) The un-renormalized NN P-wave phase shifts as a function of the laboratory kinetic energy that result 
from a long-range potential of OPE plus dimensionally regularized TPE at NNLO. The phase shifts are shown for cutoffs A 
ranging from 0.5 to 2 GeV. The Nijmegen phase-shift analysis [131 is indicated by the open diamonds. 
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FIG. 11: (Color online) The NN P-wave phase shifts as a function of the laboratory kinetic energy that result from the use 
of NNLO DR TPE and one subtraction. Here the cutoff range shown is 0.6-10 GeV. The input value of af/ is taken from 
Ref. [4l[. The Nijmegen phase-shift analysis [i^l is indicated by the open diamonds. 
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FIG. 12: (Color online) The NN P-wave phase shifts as a function of the laboratory kinetic energy that result from choosing 
v^^ = vi^ +V2Tr, with the latter chosen to be the DR TPE at NNLO, and implementing one subtraction. Here the generalized 
scattering lengths af/ are adjusted to give the best fit in the region Tiab < 100 MeV. The Nijmegen phase-shift analysis [i^l 
is indicated by the open diamonds. 
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FIG. 13: (Color online) The renormalized NN P-wave phase shifts as a function of the laboratory kinetic energy resulting when 
v^^ is chosen to be the OPE plus TPE in NNLO, where SFR is employed for the 0{Q^) part of TPE. Here the generalized 
scattering lengths aff are adjusted to give the best fit in the region Tiab < 100 MeV. The Nijmegen phase-shift analysis [i^l 
is indicated by the open diamonds. 
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FIG. 14: (Color online) The NN P-wave phase shifts as a function of the laboratory kinetic energy that result from choosing 
v^^ = WiTT + V2^, with the latter chosen to be the SFR TPE up to NNLO, and implementing one subtraction. Here the 
generalized scattering lengths af/ are adjusted to give the best fit in the region Tiab < 100 MeV. The Nijmegen phase-shift 
analysis [irj l is indicated by the open diamonds. 



